

<!DOCTYPE html>
<html lang="zh-CN" data-default-color-scheme=auto>



<head>
  <meta charset="UTF-8">
  <link rel="apple-touch-icon" sizes="76x76" href="/img/xiaotubiao.png">
  <link rel="icon" href="/img/xiaotubiao.png">
  <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=5.0, shrink-to-fit=no">
  <meta http-equiv="x-ua-compatible" content="ie=edge">
  
  <meta name="theme-color" content="#2f4154">
  <meta name="author" content="Qin Zk">
  <meta name="keywords" content="">
  
    <meta name="description" content="https:&#x2F;&#x2F;github.com&#x2F;QInzhengk&#x2F;Math-Model-and-Machine-Learning 题目及摘要关键词：数据清洗 线性插值 线性最小二乘拟合 多元回归分析 时间序列预测 一、任务重述	1微信公众号：数学建模与人工智能 二、任务的分析	2任务一的分析:任务一是一个关于数据处理的问题，要求我们对于给定原始数据，识别异常、冗余或丢失等问题数据，并完成对问题数据的修正及">
<meta property="og:type" content="article">
<meta property="og:title" content="河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测">
<meta property="og:url" content="http://example.com/2022/05/20/%E6%B2%B3%E5%8C%97%E7%9C%81%E7%AC%AC%E4%B8%89%E5%B1%8A%E7%A0%94%E7%A9%B6%E7%94%9F%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1B%E9%A2%98%EF%BC%88%E4%BA%8C%E7%AD%89%EF%BC%89%E4%BA%A4%E9%80%9A%E6%A3%80%E6%B5%8B%E5%99%A8%E6%95%B0%E6%8D%AE%E8%B4%A8%E9%87%8F%E6%8E%A7%E5%88%B6%E5%8F%8A%E9%A2%84%E6%B5%8B/index.html">
<meta property="og:site_name" content="幸运藏在努力里">
<meta property="og:description" content="https:&#x2F;&#x2F;github.com&#x2F;QInzhengk&#x2F;Math-Model-and-Machine-Learning 题目及摘要关键词：数据清洗 线性插值 线性最小二乘拟合 多元回归分析 时间序列预测 一、任务重述	1微信公众号：数学建模与人工智能 二、任务的分析	2任务一的分析:任务一是一个关于数据处理的问题，要求我们对于给定原始数据，识别异常、冗余或丢失等问题数据，并完成对问题数据的修正及">
<meta property="og:locale" content="zh_CN">
<meta property="og:image" content="https://img-blog.csdnimg.cn/20210508191330221.png">
<meta property="og:image" content="https://img-blog.csdnimg.cn/20210508191752890.png">
<meta property="og:image" content="https://img-blog.csdnimg.cn/20210508191808287.png">
<meta property="og:image" content="https://img-blog.csdnimg.cn/20210508191959475.png">
<meta property="article:published_time" content="2022-05-20T03:15:42.957Z">
<meta property="article:modified_time" content="2022-07-30T08:52:35.191Z">
<meta property="article:author" content="Qin Zk">
<meta property="article:tag" content="数学建模">
<meta name="twitter:card" content="summary_large_image">
<meta name="twitter:image" content="https://img-blog.csdnimg.cn/20210508191330221.png">
  
  
    <meta name="referrer" content="no-referrer-when-downgrade">
  
  
  <title>河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测 - 幸运藏在努力里</title>

  <link  rel="stylesheet" href="https://lib.baomitu.com/twitter-bootstrap/4.6.1/css/bootstrap.min.css" />



  <link  rel="stylesheet" href="https://lib.baomitu.com/github-markdown-css/4.0.0/github-markdown.min.css" />

  <link  rel="stylesheet" href="https://lib.baomitu.com/hint.css/2.7.0/hint.min.css" />

  <link  rel="stylesheet" href="https://lib.baomitu.com/fancybox/3.5.7/jquery.fancybox.min.css" />



<!-- 主题依赖的图标库，不要自行修改 -->
<!-- Do not modify the link that theme dependent icons -->

<link rel="stylesheet" href="//at.alicdn.com/t/font_1749284_hj8rtnfg7um.css">



<link rel="stylesheet" href="//at.alicdn.com/t/font_1736178_lbnruvf0jn.css">


<link  rel="stylesheet" href="/css/main.css" />


  <link id="highlight-css" rel="stylesheet" href="/css/highlight.css" />
  
    <link id="highlight-css-dark" rel="stylesheet" href="/css/highlight-dark.css" />
  




  <script id="fluid-configs">
    var Fluid = window.Fluid || {};
    Fluid.ctx = Object.assign({}, Fluid.ctx)
    var CONFIG = {"hostname":"example.com","root":"/","version":"1.9.0","typing":{"enable":true,"typeSpeed":70,"cursorChar":"_","loop":false,"scope":[]},"anchorjs":{"enable":true,"element":"h1,h2,h3,h4,h5,h6","placement":"left","visible":"hover","icon":""},"progressbar":{"enable":true,"height_px":3,"color":"#29d","options":{"showSpinner":false,"trickleSpeed":100}},"code_language":{"enable":true,"default":"TEXT"},"copy_btn":true,"image_caption":{"enable":true},"image_zoom":{"enable":true,"img_url_replace":["",""]},"toc":{"enable":true,"placement":"left","headingSelector":"h1,h2,h3,h4,h5,h6","collapseDepth":0},"lazyload":{"enable":true,"loading_img":"/img/loading.gif","onlypost":false,"offset_factor":2},"web_analytics":{"enable":"ture","follow_dnt":true,"baidu":null,"google":null,"gtag":null,"tencent":{"sid":null,"cid":null},"woyaola":null,"cnzz":null,"leancloud":{"app_id":null,"app_key":null,"server_url":null,"path":"window.location.pathname","ignore_local":false}},"search_path":"/local-search.xml"};

    if (CONFIG.web_analytics.follow_dnt) {
      var dntVal = navigator.doNotTrack || window.doNotTrack || navigator.msDoNotTrack;
      Fluid.ctx.dnt = dntVal && (dntVal.startsWith('1') || dntVal.startsWith('yes') || dntVal.startsWith('on'));
    }
  </script>
  <script  src="/js/utils.js" ></script>
  <script  src="/js/color-schema.js" ></script>

  
<meta name="generator" content="Hexo 6.2.0"></head>


<body>
  

  <header>
    

<div class="header-inner" style="height: 70vh;">
  <nav id="navbar" class="navbar fixed-top  navbar-expand-lg navbar-dark scrolling-navbar">
  <div class="container">
    <a class="navbar-brand" href="/">
      <strong>Qin&#39;s Blog</strong>
    </a>

    <button id="navbar-toggler-btn" class="navbar-toggler" type="button" data-toggle="collapse"
            data-target="#navbarSupportedContent"
            aria-controls="navbarSupportedContent" aria-expanded="false" aria-label="Toggle navigation">
      <div class="animated-icon"><span></span><span></span><span></span></div>
    </button>

    <!-- Collapsible content -->
    <div class="collapse navbar-collapse" id="navbarSupportedContent">
      <ul class="navbar-nav ml-auto text-center">
        
          
          
          
          
            <li class="nav-item">
              <a class="nav-link" href="/">
                <i class="iconfont icon-home-fill"></i>
                首页
              </a>
            </li>
          
        
          
          
          
          
            <li class="nav-item">
              <a class="nav-link" href="/archives/">
                <i class="iconfont icon-archive-fill"></i>
                归档
              </a>
            </li>
          
        
          
          
          
          
            <li class="nav-item">
              <a class="nav-link" href="/categories/">
                <i class="iconfont icon-category-fill"></i>
                分类
              </a>
            </li>
          
        
          
          
          
          
            <li class="nav-item">
              <a class="nav-link" href="/tags/">
                <i class="iconfont icon-tags-fill"></i>
                标签
              </a>
            </li>
          
        
          
          
          
          
            <li class="nav-item">
              <a class="nav-link" href="/about/">
                <i class="iconfont icon-user-fill"></i>
                关于
              </a>
            </li>
          
        
        
          <li class="nav-item" id="search-btn">
            <a class="nav-link" target="_self" href="javascript:;" data-toggle="modal" data-target="#modalSearch" aria-label="Search">
              &nbsp;<i class="iconfont icon-search"></i>&nbsp;
            </a>
          </li>
          
        
        
          <li class="nav-item" id="color-toggle-btn">
            <a class="nav-link" target="_self" href="javascript:;" aria-label="Color Toggle">&nbsp;<i
                class="iconfont icon-dark" id="color-toggle-icon"></i>&nbsp;</a>
          </li>
        
      </ul>
    </div>
  </div>
</nav>

  

<div id="banner" class="banner" parallax=true
     style="background: url('/img/article.png') no-repeat center center; background-size: cover;">
  <div class="full-bg-img">
    <div class="mask flex-center" style="background-color: rgba(0, 0, 0, 0.3)">
      <div class="banner-text text-center fade-in-up">
        <div class="h2">
          
            <span id="subtitle" data-typed-text="河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测"></span>
          
        </div>

        
          
  <div class="mt-3">
    
      <span class="post-meta mr-2">
        <i class="iconfont icon-author" aria-hidden="true"></i>
        Qin Zk
      </span>
    
    
      <span class="post-meta">
        <i class="iconfont icon-date-fill" aria-hidden="true"></i>
        <time datetime="2022-05-20 11:15" pubdate>
          2022年5月20日 上午
        </time>
      </span>
    
  </div>

  <div class="mt-1">
    

    

    
    
      
        <span id="busuanzi_container_page_pv" style="display: none">
          <i class="iconfont icon-eye" aria-hidden="true"></i>
          <span id="busuanzi_value_page_pv"></span> 次
        </span>
        
      
    
  </div>


        
      </div>

      
    </div>
  </div>
</div>

</div>

    <a target="_blank" rel="noopener" href="https://qzkq.github.io" class="github-corner" aria-label="View source on GitHub"><svg width="80" height="80" viewBox="0 0 250 250" style="fill:#64CEAA; color:#fff; position: absolute; top: 0; border: 0; right: 0;" aria-hidden="true"><path d="M0,0 L115,115 L130,115 L142,142 L250,250 L250,0 Z"></path><path d="M128.3,109.0 C113.8,99.7 119.0,89.6 119.0,89.6 C122.0,82.7 120.5,78.6 120.5,78.6 C119.2,72.0 123.4,76.3 123.4,76.3 C127.3,80.9 125.5,87.3 125.5,87.3 C122.9,97.6 130.6,101.9 134.4,103.2" fill="currentColor" style="transform-origin: 130px 106px;" class="octo-arm"></path><path d="M115.0,115.0 C114.9,115.1 118.7,116.5 119.8,115.4 L133.7,101.6 C136.9,99.2 139.9,98.4 142.2,98.6 C133.8,88.0 127.5,74.4 143.8,58.0 C148.5,53.4 154.0,51.2 159.7,51.0 C160.3,49.4 163.2,43.6 171.4,40.1 C171.4,40.1 176.1,42.5 178.8,56.2 C183.1,58.6 187.2,61.8 190.9,65.4 C194.5,69.0 197.7,73.2 200.1,77.6 C213.8,80.2 216.3,84.9 216.3,84.9 C212.7,93.1 206.9,96.0 205.4,96.6 C205.1,102.4 203.0,107.8 198.3,112.5 C181.9,128.9 168.3,122.5 157.7,114.1 C157.9,116.9 156.7,120.9 152.7,124.9 L141.0,136.5 C139.8,137.7 141.6,141.9 141.8,141.8 Z" fill="currentColor" class="octo-body"></path></svg></a><style>.github-corner:hover .octo-arm{animation:octocat-wave 560ms ease-in-out}@keyframes octocat-wave{0%,100%{transform:rotate(0)}20%,60%{transform:rotate(-25deg)}40%,80%{transform:rotate(10deg)}}@media (max-width:500px){.github-corner:hover .octo-arm{animation:none}.github-corner .octo-arm{animation:octocat-wave 560ms ease-in-out}}</style>
  </header>

  <main>
    
      

<div class="container-fluid nopadding-x">
  <div class="row nomargin-x">
    <div class="side-col d-none d-lg-block col-lg-2">
      
  <aside class="sidebar" style="padding-left: 2rem; margin-right: -1rem">
    <div id="toc">
  <p class="toc-header"><i class="iconfont icon-list"></i>&nbsp;目录</p>
  <div class="toc-body" id="toc-body"></div>
</div>



  </aside>


    </div>

    <div class="col-lg-8 nopadding-x-md">
      <div class="container nopadding-x-md" id="board-ctn">
        <div id="board">
          <article class="post-content mx-auto">
            <!-- SEO header -->
            <h1 style="display: none">河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测</h1>
            
            <div class="markdown-body">
              
              <p><a target="_blank" rel="noopener" href="https://github.com/QInzhengk/Math-Model-and-Machine-Learning">https://github.com/QInzhengk/Math-Model-and-Machine-Learning</a></p>
<h1 id="题目及摘要"><a href="#题目及摘要" class="headerlink" title="题目及摘要"></a><a target="_blank" rel="noopener" href="https://mp.weixin.qq.com/s?__biz=MzI5MTY1MzU1Mg==&mid=100000171&idx=1&sn=388748c28d3f18c42926be8942582872&scene=19#wechat_redirect">题目及摘要</a></h1><p><strong>关键词：数据清洗 线性插值 线性最小二乘拟合 多元回归分析 时间序列预测</strong></p>
<h1 id="一、任务重述1"><a href="#一、任务重述1" class="headerlink" title="一、任务重述	1"></a>一、任务重述	1</h1><p><a target="_blank" rel="noopener" href="https://mp.weixin.qq.com/s?__biz=MzI5MTY1MzU1Mg==&mid=100000171&idx=1&sn=388748c28d3f18c42926be8942582872&scene=19#wechat_redirect">微信公众号：数学建模与人工智能</a></p>
<h1 id="二、任务的分析2"><a href="#二、任务的分析2" class="headerlink" title="二、任务的分析	2"></a>二、任务的分析	2</h1><p>任务一的分析:任务一是一个关于数据处理的问题，要求我们对于给定原始数据，识别异常、冗余或丢失等问题数据，并完成对问题数据的修正及缺失数据的填补。对此，首先观察所给出的原始数据，发现在时间存在大量的缺失值。其次，还可能存在着冗余、异常或丢失数据。通过对原始数据的特征分析，清洗数据，对冗余值进行处理，利用标准差法和箱线图法对异常数据进行修正，接着，利用线性插值对缺失值进行填补。</p>
<p>任务二的分析：任务二可以分解成两部分，第一部分是求解两个变量间的关系，第二部分是求解多变量之间的关系，并对data 2中的缺失值进行填补。第一部分的自变量为时间，因变量为流量、速度和时间占有率三个参数。对于大样本的数据，通常采用线性拟合的办法，获得一个最佳拟合函数，表达出因变量随自变量的变化规律。第二部分的因变量可以为时间占有率，那么自变量即是流量和速度。探索三个变量之间的关系可以使用多元回归分析，得到三个变量的回归系数，进一步得到三个变量间的数学模型。最后，利用得到的数学模型对给定数据样本中的缺失值进行填补。</p>
<p>任务三的分析：任务三要求我们根据交通系统的特点建立适用、合理的预测模型，并将data 3中缺失的速度、流量和时间占有率的数据进行补充。首先，对于基于时间按照先后顺序排列的某大样本的某种统计指标的预测问题，常用的方法有神经网络预测方法和时间序列预测等方法。其中，时间序列预测模型包括平稳时间序列模型和非平稳时间序列模型。在任务二中，探讨了交通系统在时间维度上具有较强的依赖性和周期性，因此，可以合理的假设，交通系统在更短的时间内也可能具备一定的周期性。因此，针对于本题，考虑使用非平稳时间序列预测。</p>
<h1 id="三、模型假设3"><a href="#三、模型假设3" class="headerlink" title="三、模型假设	3"></a>三、模型假设	3</h1><p>（1）数据样本中不存在不精确数据；</p>
<p>（2）该路段可能有交通事故，或车辆超速行驶等情况的发生；</p>
<p>（3）该地区，公路的最高车速限制为60 km⋅h^-1；</p>
<p>（4）该地区城市道路和公路的设计通行能力为1000-4000 veh⋅h-^1；</p>
<h1 id="四、符号说明3"><a href="#四、符号说明3" class="headerlink" title="四、符号说明	3"></a>四、符号说明	3</h1><table>
<thead>
<tr>
<th>符号</th>
<th>含义</th>
</tr>
</thead>
<tbody><tr>
<td>vbef</td>
<td>代表原始数据中的速度</td>
</tr>
<tr>
<td>qbef</td>
<td>代表原始数据中的流量</td>
</tr>
<tr>
<td>v</td>
<td>代表速度</td>
</tr>
<tr>
<td>q</td>
<td>代表流量</td>
</tr>
<tr>
<td>O</td>
<td>代表占用率</td>
</tr>
<tr>
<td>outlier</td>
<td>回归系数</td>
</tr>
<tr>
<td>R</td>
<td>复相关系数</td>
</tr>
<tr>
<td>d</td>
<td>差分次数</td>
</tr>
<tr>
<td>p</td>
<td>代自由回归项数</td>
</tr>
<tr>
<td>q</td>
<td>滑动平均项数</td>
</tr>
</tbody></table>
<h1 id="五、模型的建立与求解3"><a href="#五、模型的建立与求解3" class="headerlink" title="五、模型的建立与求解	3"></a>五、模型的建立与求解	3</h1><h2 id="5-1-任务一：数据清洗与缺失值填补3"><a href="#5-1-任务一：数据清洗与缺失值填补3" class="headerlink" title="5.1 任务一：数据清洗与缺失值填补	3"></a>5.1 任务一：数据清洗与缺失值填补	3</h2><h3 id="5-1-1求解思路3"><a href="#5-1-1求解思路3" class="headerlink" title="5.1.1求解思路	3"></a>5.1.1求解思路	3</h3><p>任务一要求将一个给定的交通数据样本中存在的问题数据进行修正，并对缺失值进行填补。问题数据样本一般主要包含三种“脏数据”，分别是丢失数据、冗余数据和异常数据。首先，在统一时间步长的基础上，要建立问题数据的判别规则，分别筛选出样本中的丢失数据，冗余数据和异常数据，并相应地对其进行补齐、约简和替换，保留正确数据。接着，利用缺失值填补模型，将清洗后的数据进行缺失值填补。最终，得到具有完整数据的数据样本。解题思路示意图如下图1.1所示。<br><img src="https://img-blog.csdnimg.cn/20210508191330221.png" srcset="/img/loading.gif" lazyload alt="图1.1 数据预处理流程图"></p>
<p>图1.1 数据预处理流程图</p>
<h3 id="5-1-2数据清洗4"><a href="#5-1-2数据清洗4" class="headerlink" title="5.1.2数据清洗	4"></a>5.1.2数据清洗	4</h3><p>5.1.2.1冗余数据与丢失数据清洗模型的建立与求解</p>
<p>数据几乎没有完美的，大多数数据都存在着质量问题，因此，数据质量分析是数据分析中的至关重要的一步，是之后进行一切探索的基础。[1]对于实际生活中交通监测，在ITS数据收集过程中，由于信号不稳定及路面状况复杂等原因，得到的数据经常出现异常、丢失、冗余等问题。比如，在较少车流量的情况下，占用率却偏高，这种不符合实际情况以及不复合期望范围内的现象，便是数据异常；其次，它也违反了交通流理论。而当车辆密度过大，或者交通检测器扫描不稳定，信号传输、存储设备出现故障等情况发生时，则经常导致数据的丢失。数据冗余是指，当某段监测路段或交叉路口布置的交通检测器密度过高，导致对同一单位的重复监测，从而引起的某时间点集合中的数据具有高相似性的现象。当这些“脏数据”包含在提供到交通调控应用的数据时，必定会引起巨大的安全隐患。因此，在利用原始数据进行分析处理时，必须要先清洗数据，即，补齐丢失数据，修正异常或冗余的数据。[2]<br>首先，观察原始数据，根据每5 min取值，取时时长从1月1日00:00到1月31日00:00，共30天，因此，理论上应共有30<em>24</em>60&#x2F;5&#x3D;8640组数据，而监测到仅有8636组数据，由此可见，样本丢失了4组数据。将一段时间内收集得到数据划分为一个时间段，例如，将08:00:00到08:05:00中得到数据都定义为08：00这一时段的数据，利用python对给定数据中所有时间段进行扫描，检索是否有数据丢失，如有丢失，则进行均值填补。扫描结果得出，data 1中存在4组丢失数据，分别为：1月9日0:10，1月13日10:40，1月26日13:40和1月26日13:45。<br>接着，对原始数据进行值分析，这有利于在总体上展示数据的自然分布情况，例如，数据的唯一值、变量中的空值等。此模型中，具体统计信息如表1.1所示。</p>
<p>表格 1.1 数据统计描述</p>
<table>
<thead>
<tr>
<th>数值名称</th>
<th>速度（v）</th>
<th>流量（q）</th>
<th>占用率(O)</th>
</tr>
</thead>
<tbody><tr>
<td>总记录数</td>
<td>8640</td>
<td>8640</td>
<td>8640</td>
</tr>
<tr>
<td>非空值数</td>
<td>5981</td>
<td>6045</td>
<td>6045</td>
</tr>
<tr>
<td>空值占比</td>
<td>30.77%</td>
<td>30.03%</td>
<td>30.03%</td>
</tr>
<tr>
<td>平均值</td>
<td>65.89</td>
<td>63.68</td>
<td>10.32</td>
</tr>
<tr>
<td>标准差</td>
<td>29.96</td>
<td>39.03</td>
<td>10.00</td>
</tr>
<tr>
<td>最小值</td>
<td>3.55</td>
<td>0</td>
<td>0</td>
</tr>
<tr>
<td>下四分位数</td>
<td>52.70</td>
<td>30.00</td>
<td>3.50</td>
</tr>
<tr>
<td>中值</td>
<td>57.61</td>
<td>64.00</td>
<td>8.24</td>
</tr>
<tr>
<td>上四分位数</td>
<td>60.74</td>
<td>91.00</td>
<td>12.60</td>
</tr>
<tr>
<td>最大值</td>
<td>142.74</td>
<td>182.15</td>
<td>107.65</td>
</tr>
</tbody></table>
<p>由上表可以发现，三个参数的空值占比均为30%左右，处于5%~85%范围内，说明缺失值占据了总记录数的较大比例，如果使用删除法会导致样本量减少，削弱统计功效。因此，应当选择替换法对缺失值进行填补。此外，原始数据中可能存在着冗余数据，冗余数据的存在会影响数据分析和挖掘结果的准确性。因此，接下来对数据进行重复观测处理，在Pandas中使用duplicated方法，返回数据行每一行的检验结果，即每一行返回一个bool值，使用drop. Duplicates方法移除重复值。结果显示，给定的原始数据中不存在冗余数据。因此，接下来对原始数据进行清洗。</p>
<p>5.1.2.2异常数据清洗模型的建立与求解</p>
<p> <strong>1. 箱线图法</strong></p>
<p>关于异常值的检测与处理，箱线图是一种很常用的方法。箱线图是指主要利用一组数据的上、下边缘，中位数和上、下四分位数来展示一组数据的离群信息以及数据是否对称等信息的一种统计方法。应当注意，离群数据不一定都是异常数据。这里，箱线图法的判别公式以及判别过程如图1.2所示，其中，Q1为下四分数（25%），Q3为上四分数（75%），IQR为上四分位数与下四分位数的差。<img src="https://img-blog.csdnimg.cn/20210508191752890.png" srcset="/img/loading.gif" lazyload alt="在这里插入图片描述"><br>图1.2 箱线图判别流程图</p>
<p>由上箱线图法，结合python，对给定的原始数据进行扫描，得到结果发现速度、流量和占用率的给定原始数据都存在异常值。速度、流量和占用率的直方图和核密度图如图1.3所示，图中的直方图和核密度曲线显示，数据分布形状是有偏的。<br><img src="https://img-blog.csdnimg.cn/20210508191808287.png" srcset="/img/loading.gif" lazyload alt="在这里插入图片描述"><br>图1.3 直方图和核密度曲线</p>
<p>最后，给出异常值处理前后的统计描述对比，如表1.2所示。根据箱线图法得到速度、容量和占用率的上限临界值分别为72.80，182.5和26.25，说明此模型对异常数据的修正标定了一个合理的上限。此外，可以发现异常值的替换，不仅改变了给定的原始数据的均值，标准差和最大值，并且这些值改变后都降低了，说明此箱线图法对异常值的修正具有效果的良好。例如，给定原始数据中的最大占用率为107.65%，这明显不符合实际情况，属于异常数据。经替换后，现占用率为26.17%，符合客观事实。<br>表1.2 异常值处理前后的统计描述对比<br>    count	mean	std	min	25%	50%	75%	max	上限临界<br>speed	替换前	5981	65.89	29.96	3.55	52.70	57.61	60.74	142.74	72.80<br>    替换后	5981	55.10	11.57	3.55	52.70	57.61	60.74	67.4	<br>volume	替换前	6045	63.68	39.03	0	30.00	64.00	91.00	182.15	182.5<br>    替换后	6045	63.68	39.03	0	30.00	64.00	91.00	182.15	<br>occupancy	替换前	6045	10.21	10.00	0	3.50	8.24	12.60	107.65	26.25<br>    替换后	6045	9.45	7.28	0	3.50	8.24	12.60	26.17	</p>
<p><strong>2.阈值理论法</strong></p>
<p>此外，根据交通流理论，地点平均速度v为0，流量q不为0；流量q为0，但占用率O和地点平均速度v同时不为0；占用率为0，流量大于设定值。[3]对于交通流模型，在阈值理论中，平均速度v的取值范围为0 ≤ v ≤ fv⋅v1，其中，v1为道路的限制速度，单位为km⋅h^-1，fv为修正系数，一般取1.25<del>1.5；流量q的取值范围为0 ≤ q≤ fc⋅C⋅T&#x2F;60，其中，C为道路通行能力，单位为veh⋅h-1; T为数据采集的时间间隔,单位为min；fc为修正系数，一般取1.25</del>1.5；占用率O的取值范围为0≤ O ≤100%。[4]这里，根据城市道路一般限速为60 km⋅h^-1，取修正系数fv&#x3D;1.25，可得速度阈值为75 km⋅h^-1；道路通行能力取1520 veh⋅h^-1，fc取1.5，T为5 min，可得q的阈值为190 veh；根据实际情况推理，若路面保持畅通，则占用率O应处于一个合理范围，此处，将O阈值取为90%。根据以上得到的速度、容量和占用率阈值，设定异常值清洗规则，得到结果如图1.4和表1.3所示。由图1.4可以发现，由实际情况通过计算出发得到的直方图和核密度曲线与由箱线图法得到的结果基本相似，这说明两种方法得到的结果相似度高，模型具有较高的准确性。<br><img src="https://img-blog.csdnimg.cn/20210508191959475.png" srcset="/img/loading.gif" lazyload alt="在这里插入图片描述"><br>图1.4 直方图和核密度曲线</p>
<p>表1.3 异常值处理前后的统计描述对比</p>
<p>count		mean	std	min	25%	50%	75%	max<br>speed	替换前	5981	65.88	29.96	3.55	52.70	57.61	60.74	142.74<br>    替换后	5981	56.43	13.17	3.55	52.70	57.61	60.74	74.06<br>volume	替换前	6045	63.68	39.03	0	30.00	64.00	91.00	182.15<br>    替换后	6045	63.68	39.03	0	30.00	64.00	91.00	182.15<br>occupancy	替换前	6045	10.21	10.00	0	3.50	8.24	12.60	107.65<br>    替换后	6045	10.30	9.85	0	3.50	8.24	12.60	86.93</p>
<h3 id="5-1-3缺失值的线性插值模型的建立与求解8"><a href="#5-1-3缺失值的线性插值模型的建立与求解8" class="headerlink" title="5.1.3缺失值的线性插值模型的建立与求解	8"></a>5.1.3缺失值的线性插值模型的建立与求解	8</h3><p>经过对原始数据的清洗，数据的准确性得到了很大的提高，在此基础上，现将对缺失数据进行填补。对于缺失数据的填补，常用的方法为插值法，包括分段线性插值和分段二次插值等方法。</p>
<h2 id="5-2任务二：交通参数-时间线性拟合模型和交通参数多元回归模型9"><a href="#5-2任务二：交通参数-时间线性拟合模型和交通参数多元回归模型9" class="headerlink" title="5.2任务二：交通参数-时间线性拟合模型和交通参数多元回归模型	9"></a>5.2任务二：交通参数-时间线性拟合模型和交通参数多元回归模型	9</h2><h3 id="5-2-1求解思路9"><a href="#5-2-1求解思路9" class="headerlink" title="5.2.1求解思路	9"></a>5.2.1求解思路	9</h3><h3 id="5-2-2数据清洗模型的建立与求解9"><a href="#5-2-2数据清洗模型的建立与求解9" class="headerlink" title="5.2.2数据清洗模型的建立与求解	9"></a>5.2.2数据清洗模型的建立与求解	9</h3><h3 id="5-2-3交通参数-时间线性最小二乘拟合模型的建立与求解10"><a href="#5-2-3交通参数-时间线性最小二乘拟合模型的建立与求解10" class="headerlink" title="5.2.3交通参数-时间线性最小二乘拟合模型的建立与求解	10"></a>5.2.3交通参数-时间线性最小二乘拟合模型的建立与求解	10</h3><h3 id="5-2-4多元回归分析模型的建立与求解16"><a href="#5-2-4多元回归分析模型的建立与求解16" class="headerlink" title="5.2.4多元回归分析模型的建立与求解	16"></a>5.2.4多元回归分析模型的建立与求解	16</h3><h3 id="5-2-5数据填补21"><a href="#5-2-5数据填补21" class="headerlink" title="5.2.5数据填补	21"></a>5.2.5数据填补	21</h3><h2 id="5-3-任务三：非平稳时间序列预测模型22"><a href="#5-3-任务三：非平稳时间序列预测模型22" class="headerlink" title="5.3 任务三：非平稳时间序列预测模型	22"></a>5.3 任务三：非平稳时间序列预测模型	22</h2><h3 id="5-3-1求解思路22"><a href="#5-3-1求解思路22" class="headerlink" title="5.3.1求解思路	22"></a>5.3.1求解思路	22</h3><h3 id="5-3-2-ARIMA模型的建立与求解22"><a href="#5-3-2-ARIMA模型的建立与求解22" class="headerlink" title="5.3.2 ARIMA模型的建立与求解	22"></a>5.3.2 ARIMA模型的建立与求解	22</h3><h1 id="六、模型评价与改进26"><a href="#六、模型评价与改进26" class="headerlink" title="六、模型评价与改进	26"></a>六、模型评价与改进	26</h1><h1 id="七、参考文献27"><a href="#七、参考文献27" class="headerlink" title="七、参考文献	27"></a>七、参考文献	27</h1><h1 id="八、附录28"><a href="#八、附录28" class="headerlink" title="八、附录	28"></a>八、附录	28</h1><figure class="highlight python"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br><span class="line">45</span><br><span class="line">46</span><br><span class="line">47</span><br><span class="line">48</span><br><span class="line">49</span><br><span class="line">50</span><br><span class="line">51</span><br><span class="line">52</span><br><span class="line">53</span><br><span class="line">54</span><br><span class="line">55</span><br><span class="line">56</span><br><span class="line">57</span><br><span class="line">58</span><br><span class="line">59</span><br><span class="line">60</span><br><span class="line">61</span><br><span class="line">62</span><br><span class="line">63</span><br><span class="line">64</span><br><span class="line">65</span><br><span class="line">66</span><br><span class="line">67</span><br><span class="line">68</span><br><span class="line">69</span><br><span class="line">70</span><br><span class="line">71</span><br><span class="line">72</span><br><span class="line">73</span><br><span class="line">74</span><br><span class="line">75</span><br><span class="line">76</span><br><span class="line">77</span><br><span class="line">78</span><br><span class="line">79</span><br><span class="line">80</span><br><span class="line">81</span><br><span class="line">82</span><br><span class="line">83</span><br><span class="line">84</span><br><span class="line">85</span><br><span class="line">86</span><br><span class="line">87</span><br><span class="line">88</span><br><span class="line">89</span><br><span class="line">90</span><br><span class="line">91</span><br><span class="line">92</span><br><span class="line">93</span><br><span class="line">94</span><br><span class="line">95</span><br><span class="line">96</span><br><span class="line">97</span><br><span class="line">98</span><br></pre></td><td class="code"><pre><code class="hljs python"><span class="hljs-keyword">import</span> pandas <span class="hljs-keyword">as</span> pd<br>a=pd.read_excel(<span class="hljs-string">&quot;end12.xlsx&quot;</span>,usecols=<span class="hljs-built_in">range</span>(<span class="hljs-number">1</span>,<span class="hljs-number">4</span>))  <span class="hljs-comment">#提取第2列到第4列的数据</span><br>b=a.values  <span class="hljs-comment">#提取其中的数据</span><br>c=a.describe()  <span class="hljs-comment">#对数据进行统计描述</span><br><span class="hljs-built_in">print</span>(c)<br><br><br><span class="hljs-keyword">import</span> pandas <span class="hljs-keyword">as</span> pd<br>a=pd.read_excel(<span class="hljs-string">&quot;data 1.xlsx&quot;</span>)<br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;是否存在重复观测：&quot;</span>,<span class="hljs-built_in">any</span>(a.duplicated()))  <span class="hljs-comment">#输出：True</span><br>a.drop_duplicates(inplace=<span class="hljs-literal">True</span>)  <span class="hljs-comment">#inplace=True时，直接删除a中的重复数据</span><br>f=pd.ExcelWriter(<span class="hljs-string">&#x27;data_1.xlsx&#x27;</span>)  <span class="hljs-comment">#创建文件对象</span><br>a.to_excel(f)  <span class="hljs-comment">#把a写入新Excel文件中</span><br>f.save()       <span class="hljs-comment">#保存文件，数据才真正写入Excel文件</span><br><br><br><span class="hljs-keyword">from</span> pandas <span class="hljs-keyword">import</span> read_excel<br><span class="hljs-keyword">import</span> pandas <span class="hljs-keyword">as</span> pd<br>a=read_excel(<span class="hljs-string">&quot;data 2.xlsx&quot;</span>)<br>b=a.fillna(value=&#123;<span class="hljs-string">&#x27;speed&#x27;</span>:a.speed.interpolate(),<br>                   <span class="hljs-string">&#x27;volume&#x27;</span>:a.volume.interpolate(),<br>                   <span class="hljs-string">&#x27;occupancy&#x27;</span>:a.occupancy.interpolate()&#125;) <span class="hljs-comment">#使用线性插值替换</span><br>f=pd.ExcelWriter(<span class="hljs-string">&#x27;data22.xlsx&#x27;</span>)  <span class="hljs-comment">#创建文件对象</span><br>b.to_excel(f)  <span class="hljs-comment">#把a写入新Excel文件中</span><br>f.save()       <span class="hljs-comment">#保存文件，数据才真正写入Excel文件</span><br><br><br><span class="hljs-keyword">from</span> pandas <span class="hljs-keyword">import</span> read_csv<br><span class="hljs-keyword">import</span> numpy <span class="hljs-keyword">as</span> np<br><span class="hljs-keyword">import</span> pandas <span class="hljs-keyword">as</span> pd<br><span class="hljs-keyword">import</span> matplotlib.pyplot <span class="hljs-keyword">as</span> plt<br>a=read_csv(<span class="hljs-string">&quot;data 1.csv&quot;</span>)<br>mu=a.occupancy.mean()<br>s=a.occupancy.std()<br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;标准差法异常值上限检测：&quot;</span>,<span class="hljs-built_in">any</span>(a.occupancy&gt;mu+<span class="hljs-number">2</span>*s))  <span class="hljs-comment">#输出：True</span><br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;标准差法异常值下限检测：&quot;</span>,<span class="hljs-built_in">any</span>(a.occupancy&lt;mu-<span class="hljs-number">2</span>*s))  <span class="hljs-comment">#输出：False</span><br>Q1=a.occupancy.quantile(<span class="hljs-number">0.25</span>)  <span class="hljs-comment">#计算下四分位数</span><br>Q3=a.occupancy.quantile(<span class="hljs-number">0.75</span>)  <span class="hljs-comment">#计算上四分位数</span><br>IQR=Q3-Q1<br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;箱线图法异常值上限检测：&quot;</span>,<span class="hljs-built_in">any</span>(a.occupancy&gt;Q3+<span class="hljs-number">1.5</span>*IQR))  <span class="hljs-comment">#输出：True</span><br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;箱线图法异常值下限检测：&quot;</span>,<span class="hljs-built_in">any</span>(a.occupancy&lt;Q1-<span class="hljs-number">1.5</span>*IQR))  <span class="hljs-comment">#输出：False</span><br>plt.style.use(<span class="hljs-string">&#x27;ggplot&#x27;</span>)  <span class="hljs-comment">#设置绘图风格</span><br>a.occupancy.plot(kind=<span class="hljs-string">&#x27;hist&#x27;</span>,bins=<span class="hljs-number">30</span>,density=<span class="hljs-literal">True</span>)  <span class="hljs-comment">#绘制直方图</span><br>a.occupancy.plot(kind=<span class="hljs-string">&#x27;kde&#x27;</span>)  <span class="hljs-comment">#绘制核密度曲线</span><br>plt.show()<br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;异常值替换前的数据统计特征&quot;</span>,a.occupancy.describe())<br>UB=Q3+<span class="hljs-number">1.5</span>*IQR;<br>st=a.occupancy[a.occupancy&lt;UB].<span class="hljs-built_in">max</span>()  <span class="hljs-comment">#找出低于判别上限的最大值</span><br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;判别异常值的上限临界值为:&quot;</span>,UB)<br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;用以替换异常值的数据为：&quot;</span>,st)<br>a.loc[a.occupancy&gt;UB, <span class="hljs-string">&#x27;occupancy&#x27;</span>]=st  <span class="hljs-comment">#替换超过判别上限异常值</span><br><span class="hljs-built_in">print</span>(<span class="hljs-string">&quot;异常值替换后的数据统计特征&quot;</span>,a.occupancy.describe())<br>f=pd.ExcelWriter(<span class="hljs-string">&#x27;data_11.xlsx&#x27;</span>)  <span class="hljs-comment">#创建文件对象</span><br>a.to_excel(f)  <span class="hljs-comment">#把a写入新Excel文件中</span><br>f.save()       <span class="hljs-comment">#保存文件，数据才真正写入Excel文件</span><br><br><br>%%时间转换成自变量<br>clc<br>clear<br>close <span class="hljs-built_in">all</span><br>%%时间间隔转换<br>% DAY=xlsread(<span class="hljs-string">&#x27;data 3.csv&#x27;</span>,<span class="hljs-string">&#x27;data 3&#x27;</span>,<span class="hljs-string">&#x27;I1:I8588&#x27;</span>);<br>% HOUR=xlsread(<span class="hljs-string">&#x27;data 3.csv&#x27;</span>,<span class="hljs-string">&#x27;data 3&#x27;</span>,<span class="hljs-string">&#x27;J1:J8588&#x27;</span>);<br>% MINTURE=xlsread(<span class="hljs-string">&#x27;data 3.csv&#x27;</span>,<span class="hljs-string">&#x27;data 3&#x27;</span>,<span class="hljs-string">&#x27;K1:K8588&#x27;</span>);<br>% k=<span class="hljs-number">1</span>;<br>% <span class="hljs-keyword">for</span> i=<span class="hljs-number">2</span>:<span class="hljs-number">1</span>:<span class="hljs-number">8588</span><br>%     D=DAY(i,<span class="hljs-number">1</span>);D_1=DAY(i-<span class="hljs-number">1</span>,<span class="hljs-number">1</span>);<br>%     H=HOUR(i,<span class="hljs-number">1</span>);<br>%     M=MINTURE(i,<span class="hljs-number">1</span>);<br>%     <span class="hljs-keyword">if</span> D~=D_1<br>%         k=k+<span class="hljs-number">1</span>;  <br>%     end<br>%     X(<span class="hljs-number">1</span>,<span class="hljs-number">1</span>)=MINTURE(<span class="hljs-number">1</span>,<span class="hljs-number">1</span>); X(<span class="hljs-number">1</span>,<span class="hljs-number">2</span>)=<span class="hljs-number">1</span>;<br>%     X(i,<span class="hljs-number">2</span>)=k;<br>%     X(i,<span class="hljs-number">1</span>)=<span class="hljs-number">60</span>*H+M;<br>% end<br>% <span class="hljs-keyword">for</span> i=<span class="hljs-number">1</span>:<span class="hljs-number">1</span>:<span class="hljs-number">8588</span><br>%     X(i,<span class="hljs-number">1</span>)=X(i,<span class="hljs-number">1</span>)/<span class="hljs-number">5</span>;<br>% end<br>% %曲线拟合<br>% 读取数据<br>time=xlsread(<span class="hljs-string">&#x27;time1.xlsx&#x27;</span>,<span class="hljs-string">&#x27;time1&#x27;</span>,<span class="hljs-string">&#x27;A2:A5188&#x27;</span>);<br>speed=xlsread(<span class="hljs-string">&#x27;time1.xlsx&#x27;</span>,<span class="hljs-string">&#x27;time1&#x27;</span>,<span class="hljs-string">&#x27;B2:B5188&#x27;</span>);<br>volume=xlsread(<span class="hljs-string">&#x27;time1.xlsx&#x27;</span>,<span class="hljs-string">&#x27;time1&#x27;</span>,<span class="hljs-string">&#x27;C2:C5188&#x27;</span>);<br>occupancy=xlsread(<span class="hljs-string">&#x27;time1.xlsx&#x27;</span>,<span class="hljs-string">&#x27;time1&#x27;</span>,<span class="hljs-string">&#x27;D2:D5188&#x27;</span>);<br>p=xlsread(<span class="hljs-string">&#x27;time1.xlsx&#x27;</span>,<span class="hljs-string">&#x27;time1&#x27;</span>,<span class="hljs-string">&#x27;E2:E5188&#x27;</span>);<br>scatter3(volume,speed,occupancy,<span class="hljs-string">&#x27;r.&#x27;</span>);<br>title(<span class="hljs-string">&#x27;流量速度占有率散点图&#x27;</span>)<br>xlabel(<span class="hljs-string">&#x27;volume&#x27;</span>);<br>ylabel(<span class="hljs-string">&#x27;speed&#x27;</span>);<br>zlabel(<span class="hljs-string">&#x27;occupancy&#x27;</span>);<br>hold on<br>% scatter(time,volume,<span class="hljs-string">&#x27;r.&#x27;</span>);<br>% hold on<br>% scatter(time,occupancy,<span class="hljs-string">&#x27;g.&#x27;</span>);<br>% hold on<br>Cftool<br></code></pre></td></tr></table></figure>

<h1 id="未完待续-更新时间：2022-x2F-4-x2F-13）"><a href="#未完待续-更新时间：2022-x2F-4-x2F-13）" class="headerlink" title="未完待续(更新时间：2022&#x2F;4&#x2F;13）"></a>未完待续(更新时间：2022&#x2F;4&#x2F;13）</h1><p>资源：<a target="_blank" rel="noopener" href="https://download.csdn.net/download/qq_45832050/85180157">河北省研究生数学建模竞赛试题（第二、三、四届）</a></p>

              
            </div>
            <hr/>
            <div>
              <div class="post-metas my-3">
  
    <div class="post-meta mr-3 d-flex align-items-center">
      <i class="iconfont icon-category"></i>
      

<span class="category-chains">
  
  
    
      <span class="category-chain">
        
  <a href="/categories/%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1/" class="category-chain-item">数学建模</a>
  
  

      </span>
    
  
</span>

    </div>
  
  
    <div class="post-meta">
      <i class="iconfont icon-tags"></i>
      
        <a href="/tags/%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1/">#数学建模</a>
      
    </div>
  
</div>


              
  

  <div class="license-box my-3">
    <div class="license-title">
      <div>河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测</div>
      <div>http://example.com/2022/05/20/河北省第三届研究生数学建模B题（二等）交通检测器数据质量控制及预测/</div>
    </div>
    <div class="license-meta">
      
        <div class="license-meta-item">
          <div>作者</div>
          <div>Qin Zk</div>
        </div>
      
      
        <div class="license-meta-item license-meta-date">
          <div>发布于</div>
          <div>2022年5月20日</div>
        </div>
      
      
      <div class="license-meta-item">
        <div>许可协议</div>
        <div>
          
            
            
              <a target="_blank" href="https://creativecommons.org/licenses/by/4.0/">
              <span class="hint--top hint--rounded" aria-label="BY - 署名">
                <i class="iconfont icon-by"></i>
              </span>
              </a>
            
          
        </div>
      </div>
    </div>
    <div class="license-icon iconfont"></div>
  </div>



              
                <div class="post-prevnext my-3">
                  <article class="post-prev col-6">
                    
                    
                      <a href="/2022/07/30/ACM%E6%A8%A1%E5%BC%8F%E4%B9%8B%E8%BE%93%E5%85%A5%E8%BE%93%E5%87%BA%EF%BC%88Java_Python%E4%BE%8B%E9%A2%98%EF%BC%89/" title="ACM模式之输入输出（Java/Python例题）">
                        <i class="iconfont icon-arrowleft"></i>
                        <span class="hidden-mobile">ACM模式之输入输出（Java/Python例题）</span>
                        <span class="visible-mobile">上一篇</span>
                      </a>
                    
                  </article>
                  <article class="post-next col-6">
                    
                    
                      <a href="/2022/05/20/%E5%9F%BA%E4%BA%8EHexo%E7%9A%84%E4%B8%BB%E9%A2%98Fluid%E6%90%AD%E5%BB%BAGithub%E5%8D%9A%E5%AE%A2/" title="基于Hexo的主题Fluid搭建Github博客">
                        <span class="hidden-mobile">基于Hexo的主题Fluid搭建Github博客</span>
                        <span class="visible-mobile">下一篇</span>
                        <i class="iconfont icon-arrowright"></i>
                      </a>
                    
                  </article>
                </div>
              
            </div>

            
  <article id="comments" lazyload>
    
  <div id="lv-container" data-id="city" data-uid="MTAyMC81NjkxOC8zMzM4Mg==">
    <script type="text/javascript">
      Fluid.utils.loadComments('#lv-container', function() {
        Fluid.utils.createScript('https://cdn-city.livere.com/js/embed.dist.js');
      });
    </script>
    <noscript>Please enable JavaScript to view the comments</noscript>
  </div>


  </article>


          </article>
        </div>
      </div>
    </div>

    <div class="side-col d-none d-lg-block col-lg-2">
      

    </div>
  </div>
</div>





  



  



  



  



  






    

    
      <a id="scroll-top-button" aria-label="TOP" href="#" role="button">
        <i class="iconfont icon-arrowup" aria-hidden="true"></i>
      </a>
    

    
      <div class="modal fade" id="modalSearch" tabindex="-1" role="dialog" aria-labelledby="ModalLabel"
     aria-hidden="true">
  <div class="modal-dialog modal-dialog-scrollable modal-lg" role="document">
    <div class="modal-content">
      <div class="modal-header text-center">
        <h4 class="modal-title w-100 font-weight-bold">搜索</h4>
        <button type="button" id="local-search-close" class="close" data-dismiss="modal" aria-label="Close">
          <span aria-hidden="true">&times;</span>
        </button>
      </div>
      <div class="modal-body mx-3">
        <div class="md-form mb-5">
          <input type="text" id="local-search-input" class="form-control validate">
          <label data-error="x" data-success="v" for="local-search-input">关键词</label>
        </div>
        <div class="list-group" id="local-search-result"></div>
      </div>
    </div>
  </div>
</div>

    

    
  </main>

  <footer>
    <div class="footer-inner">
  
    <div class="footer-content">
       <a href="https://hexo.io" target="_blank" rel="nofollow noopener"><span>Hexo</span></a> <i class="iconfont icon-love"></i> <a href="https://github.com/fluid-dev/hexo-theme-fluid" target="_blank" rel="nofollow noopener"><span>Fluid</span></a> 
    </div>
  
  
    <div class="statistics">
  
  

  
    
      <span id="busuanzi_container_site_pv" style="display: none">
        总访问量 
        <span id="busuanzi_value_site_pv"></span>
         次
      </span>
    
    
      <span id="busuanzi_container_site_uv" style="display: none">
        总访客数 
        <span id="busuanzi_value_site_uv"></span>
         人
      </span>
    
    
  
</div>

  
  
  
</div>

  </footer>

  <!-- Scripts -->
  
  <script  src="https://lib.baomitu.com/nprogress/0.2.0/nprogress.min.js" ></script>
  <link  rel="stylesheet" href="https://lib.baomitu.com/nprogress/0.2.0/nprogress.min.css" />

  <script>
    NProgress.configure({"showSpinner":false,"trickleSpeed":100})
    NProgress.start()
    window.addEventListener('load', function() {
      NProgress.done();
    })
  </script>


<script  src="https://lib.baomitu.com/jquery/3.6.0/jquery.min.js" ></script>
<script  src="https://lib.baomitu.com/twitter-bootstrap/4.6.1/js/bootstrap.min.js" ></script>
<script  src="/js/events.js" ></script>
<script  src="/js/plugins.js" ></script>


  <script  src="https://lib.baomitu.com/typed.js/2.0.12/typed.min.js" ></script>
  <script>
    (function (window, document) {
      var typing = Fluid.plugins.typing;
      var subtitle = document.getElementById('subtitle');
      if (!subtitle || !typing) {
        return;
      }
      var text = subtitle.getAttribute('data-typed-text');
      
        typing(text);
      
    })(window, document);
  </script>




  

  

  

  

  

  

  




  
    <script  src="/js/img-lazyload.js" ></script>
  




  
<script>
  Fluid.utils.createScript('https://lib.baomitu.com/tocbot/4.18.0/tocbot.min.js', function() {
    var toc = jQuery('#toc');
    if (toc.length === 0 || !window.tocbot) { return; }
    var boardCtn = jQuery('#board-ctn');
    var boardTop = boardCtn.offset().top;

    window.tocbot.init({
      tocSelector     : '#toc-body',
      contentSelector : '.markdown-body',
      headingSelector : CONFIG.toc.headingSelector || 'h1,h2,h3,h4,h5,h6',
      linkClass       : 'tocbot-link',
      activeLinkClass : 'tocbot-active-link',
      listClass       : 'tocbot-list',
      isCollapsedClass: 'tocbot-is-collapsed',
      collapsibleClass: 'tocbot-is-collapsible',
      collapseDepth   : CONFIG.toc.collapseDepth || 0,
      scrollSmooth    : true,
      headingsOffset  : -boardTop
    });
    if (toc.find('.toc-list-item').length > 0) {
      toc.css('visibility', 'visible');
    }
  });
</script>


  <script>
  (function() {
    var enableLang = CONFIG.code_language.enable && CONFIG.code_language.default;
    var enableCopy = CONFIG.copy_btn;
    if (!enableLang && !enableCopy) {
      return;
    }

    function getBgClass(ele) {
      return Fluid.utils.getBackgroundLightness(ele) >= 0 ? 'code-widget-light' : 'code-widget-dark';
    }

    var copyTmpl = '';
    copyTmpl += '<div class="code-widget">';
    copyTmpl += 'LANG';
    copyTmpl += '</div>';
    jQuery('.markdown-body pre').each(function() {
      var $pre = jQuery(this);
      if ($pre.find('code.mermaid').length > 0) {
        return;
      }
      if ($pre.find('span.line').length > 0) {
        return;
      }

      var lang = '';

      if (enableLang) {
        lang = CONFIG.code_language.default;
        if ($pre[0].children.length > 0 && $pre[0].children[0].classList.length >= 2 && $pre.children().hasClass('hljs')) {
          lang = $pre[0].children[0].classList[1];
        } else if ($pre[0].getAttribute('data-language')) {
          lang = $pre[0].getAttribute('data-language');
        } else if ($pre.parent().hasClass('sourceCode') && $pre[0].children.length > 0 && $pre[0].children[0].classList.length >= 2) {
          lang = $pre[0].children[0].classList[1];
          $pre.parent().addClass('code-wrapper');
        } else if ($pre.parent().hasClass('markdown-body') && $pre[0].classList.length === 0) {
          $pre.wrap('<div class="code-wrapper"></div>');
        }
        lang = lang.toUpperCase().replace('NONE', CONFIG.code_language.default);
      }
      $pre.append(copyTmpl.replace('LANG', lang).replace('code-widget">',
        getBgClass($pre[0]) + (enableCopy ? ' code-widget copy-btn" data-clipboard-snippet><i class="iconfont icon-copy"></i>' : ' code-widget">')));

      if (enableCopy) {
        Fluid.utils.createScript('https://lib.baomitu.com/clipboard.js/2.0.10/clipboard.min.js', function() {
          var clipboard = new window.ClipboardJS('.copy-btn', {
            target: function(trigger) {
              var nodes = trigger.parentNode.childNodes;
              for (var i = 0; i < nodes.length; i++) {
                if (nodes[i].tagName === 'CODE') {
                  return nodes[i];
                }
              }
            }
          });
          clipboard.on('success', function(e) {
            e.clearSelection();
            e.trigger.innerHTML = e.trigger.innerHTML.replace('icon-copy', 'icon-success');
            setTimeout(function() {
              e.trigger.innerHTML = e.trigger.innerHTML.replace('icon-success', 'icon-copy');
            }, 2000);
          });
        });
      }
    });
  })();
</script>


  
<script>
  Fluid.utils.createScript('https://lib.baomitu.com/anchor-js/4.3.1/anchor.min.js', function() {
    window.anchors.options = {
      placement: CONFIG.anchorjs.placement,
      visible  : CONFIG.anchorjs.visible
    };
    if (CONFIG.anchorjs.icon) {
      window.anchors.options.icon = CONFIG.anchorjs.icon;
    }
    var el = (CONFIG.anchorjs.element || 'h1,h2,h3,h4,h5,h6').split(',');
    var res = [];
    for (var item of el) {
      res.push('.markdown-body > ' + item.trim());
    }
    if (CONFIG.anchorjs.placement === 'left') {
      window.anchors.options.class = 'anchorjs-link-left';
    }
    window.anchors.add(res.join(', '));
  });
</script>


  
<script>
  Fluid.utils.createScript('https://lib.baomitu.com/fancybox/3.5.7/jquery.fancybox.min.js', function() {
    Fluid.plugins.fancyBox();
  });
</script>


  <script>Fluid.plugins.imageCaption();</script>

  <script  src="/js/local-search.js" ></script>

  <script defer src="https://busuanzi.ibruce.info/busuanzi/2.3/busuanzi.pure.mini.js" ></script>





<!-- 主题的启动项，将它保持在最底部 -->
<!-- the boot of the theme, keep it at the bottom -->
<script  src="/js/boot.js" ></script>


  

  <noscript>
    <div class="noscript-warning">博客在允许 JavaScript 运行的环境下浏览效果更佳</div>
  </noscript>
</body>
</html>
